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We present a method for solving trapped few-body problems and apply it to three equal-mass 
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simplifies the calculation of energy eigenstates of the full Hamiltonian in a truncated Hilbert space 
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constructed from the trap Hamiltonian eigenstates. Particle superselection rules are implemented 
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by choosing the relevant representations of C6i;. We find that the one-dimensional system shows 
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Q\ and reducibility in this system and in standard few-body approximation techniques. 
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I. INTRODUCTION 



The advent of precision control and measurement of interacting ultracold atomic systems 
has provided a productive testing ground for few-body physics. Predictions from condensed 
matter and nuclear physics made in the 1960's and 1970's, like fermionization of strongly- 
interacting bosons and universality in Efimovian states, have been verified using ultracold 
atoms as a kind of quantum simulator. The signatures of few body physics can also be seen in 
the statistical and dynamical properties of ultracold quantum gases. Further, deterministic 
loading of a few cold atoms into deep optical lattices has made feasible the prospect of 
spectroscopic measurements of highly-correlated atomic states. 

This last application in particular has generated a lot of interest in analytical results 
and theoretical methods for a few particles that are trapped in an approximately harmonic 
well and interacting via short-range potentials p~H23]. See [23] for a recent review and a 
more complete bibliography, including experimental progress. One motivation for this line 
of research is, like untrapped few-body systems, trapped systems offer the possibility for 
observing universal properties P_2] [2^], i.e. few-body observables whose value is independent 
of the short-range interaction details or particle structure |26j . 




For harmonic traps, exact results exist for two bodies in the limit of a zero-range, con- 
tact interaction [lj [2], sometimes called the scaling limit. For three harmonically-trapped 
particles, there are analytic solutions in the joint case of the scaling limit combined with 
the unitary (or resonance) limit (H El [32], i-e strong interactions, usually described by a 
large two-body scattering length. However, except for the unitary limit, solving for the en- 
ergies and eigenstates of few-body quantum systems requires approximations. This article 
describes a method to analyze the case of three equal-mass particles in one-dimension in- 
teracting via a zero-range contact interaction. We will calculate the variation of the energy 
spectrum as a function of interaction strength as it varies between the repulsive unitary 
limit to the attractive unitary limit. We will combine techniques that have have been pre- 
viously employed with a geometrical approach to implementing the discrete symmetries of 
the Hamiltonian and any superselection rules due to indistinguishable particles. In partic- 
ular, we find that representations of the symmetry groups C6 V and C2 V allow us to present 
a unified classification scheme for the atom, dimer and trimer states with different three- 
particle content, including two and three identical bosons or fermions, as well as no identical 
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particles. This is a generalizable method that maximizes separability and reducibility in 
trapped few-body systems and therefore enhances computational efficiency Our method 
is most similar in spirit to the techniques of |27J, who use discrete symmetries to classify 
three-particle states in two dimensions, and [2T], who use discrete symmetries to reduce the 
problem of two-particles in an anisotropic well. 

A. Separability 

Approximations are particularly useful when they exploit symmetries to map the few- 
body system in a controlled way onto a related mathematical structure that has more sep- 
arability and reducibility. The standard example of separability is found in the transfor- 
mation from the particle coordinates to the center-of-mass/relative coordinates, e.g. Jacobi 
coordinates, allowing the energy eigenstates to be expressed as an unentangled product of 
center-of-mass and relative wave functions. This separation is exact for any configuration of 
free particles due to Galilean invariance, and it is also exact for a quadratic trap (even with 
unequal masses) as long as the non-interacting trap oscillation frequencies for the particles 
are all the same [9] . Here we will restrict ourselves to the case of equal masses for simplicity. 
This will ensure center-of-mass vs. relative separability even when we add interparticle in- 
teractions to construct the full Hamiltonian. The case of equal masses also provides another 
kind of separability. As we discuss below, the trap Hamiltonian has SU(3) symmetry, and 
this symmetry allows the trap Hamiltonian to be superintegrable, and therefore separable 
in multiple coordinate systems. 

Beyond these exact separability results, a standard class of approximations to increase 
separability are called adiabatic approximations. For example, the Born-Oppenheimer ap- 
proximation effectively separates the electronic and nuclear wave functions by assuming the 
electronic wave function varies at a much faster scale than the nuclear wave function. In the 
context of the three-body problem, a similar approximation is made when using hyperspher- 
ical coordinates to represent the relative degrees of freedom when solving time-independent 
Schrodinger equation J3J El US] or Faddeev equations [U |2H] in coordinate space. A good 
low-energy approximation is to effectively separate the wave function in terms of the hy- 
perradius and hyperangle by neglecting the variation of the wave function in hyperradius 
when solving the hyperangular equation. This adiabatic approximation becomes exact in 
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the unitary limit [3]. This occurs because the zero- range contact interaction preserves a 
dynamical scale invariance that introduces a hidden SO(2, 1) symmetry to the problem [§]. 

An advantage of working in one-dimension is that the six-dimensional 'three-body hy- 
perspherical relative coordinates' are just two-dimensional cylindrical coordinates. In one- 
dimension, the separation between hyperangle and hyperradius that emerges in the unitary 
limit will be easily evident; the three-fermion states are the exact solution in this limit and 
we will see that they are separable in this sense. 

B. Reducibility 

The truncated Hilbert space approximation maps the infinite-dimensional Hilbert space 
of the trapped three-body system onto a truncated model space constructed from exact 
solutions of the trap Hamiltonian without particle interactions. In a sense, this employs the 
most basic sense of reducibility: all finite-dimensional Hilbert spaces can be reduced to a sum 
of subspaces corresponding to degenerate eigenvalues of the full Hamiltonian. Procedures 
can be used to generate an effective Hamiltonian that is corrected for the energy cut-off due 
to truncation and which handle the difficulties in implementing a possibly non-Hermitian, 
non-renormalizaible, highly singular zero-range interparticle interaction in a consistent way. 
For example, such a method has employed by the no-core shell model with a Hamiltonian 
renormalized using effective field theory [T2l [T71 1291 130] . Similar truncated Hilbert space 
models but with different approaches to regularization are employed in [10] and [15j ITS] . 
In one-dimension, we will not have these difficulties with renormalization and regulariza- 
tion, but we will be able to clearly see the effect of the energy cut-off on the deviation of 
approximate results from the exact values in the universal, unitary limit. 

To further enhance reducibility, we will exploit the discrete point symmetry C$ v of the 
relative configuration space. This symmetry group is the geometrical realization of the 
permutation group of three particles, combined with the overall parity inversion symmetry. 
Utilizing properties of the irreducible representations of this group make two calculational 
steps much more efficient, as is well-known in applications from molecular and chemical 
physics. First, diagonalizing the full Hamiltonian in the truncated Hilbert space can be 
made more efficient because there are only non-zero matrix elements of the interaction 
between states in the same representation of Cq v . Second, implementing superselection rules 
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to handle indistinguishable particles is automatic: only a few representations are relevant 
depending on which particles (if any) are indistinguishable. We will apply our method to 
the cases of three indistinguishable bosons (BBB), two bosons and another particle (BBX), 
three indistinguishable fermions (FFF), two indistinguishable fermions and one other particle 
(FFX) and three distinguishable particles (XYZ). As discussed below, we assume there is no 
entanglement between the internal states of the particles and the motional state, or among 
the internal states of the particle. It is specifically the application of the Cq v symmetry, 
when used in combination with hyperspherical coordinates and a truncated Hilbert space 
approach, that we argue provides a meaningful efficiency improvement over existing methods. 
In particular, several authors have discussed the computational difficulties of implementing 
particle symmetries in a variety of truncated Hilbert space models pT5HTTl 129] , and we believe 
this is an elegant, extensible method for solving these difficulties efficiently. 

C. Motivation 

Besides our motivation to demonstrate our method and its connection to previous meth- 
ods, there is genuine theoretical and experimental interest in quantum gases and other ul- 
tracold atomic systems trapped in one-dimension. See, for example, the recent review [3"T] . 
Therefore, we do not consider this application to be merely a "toy model". An effective one- 
dimensional description emerges as a limit of the three-dimensional the harmonic 
trap becomes long in one direction and tight in the transverse plane, see for example [32J. 
A surprising number of the structures of interest that are found in the three-dimensional 
case, universal and otherwise, will also be manifest in the one-dimensional problem [33]. For 
example, when the interaction is attractive and sufficiently strong for three indistinguishable 
bosons, the system manifests all the classes of bound states found for attractive contact in- 
teractions in a three-dimensional harmonic potential: trimer states, dimer plus single atom 
states, and states of three bound atoms. In the spectra below, the emergence of universality 
in the unitary limit also will be evident. 

In the next section, we describe the model Hamiltonian for three equal-mass particles in 
one-dimension, establish notation, and begin a discussion of the symmetries of the Hamil- 
tonian that will continue into Section III. In both sections, the usefulness of the Jacobi 
hyperspherical (i.e. cylindrical) relative coordinates for visualizing and analyzing configura- 
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tion space will be emphasized. In Section IV, we provide results for the energy eigenstates 
and show how the richness of three-body physics is evidenced even in this simple model. 
Beyond the advantages of analytic power and explanatory appeal, studying this system in 
one-dimension also suggests extensions of this method to more complicated cases, like more 
particles, non-equal mass particles, genuine multi-body interactions, higher dimensions, and 
anisotropic traps. The final section contains a few remarks about these possible applications. 



II. THE MODEL 



The total Hamiltonian H can be written as the sum of the harmonic trap Hamiltonian 
H , which in position space is 

H ° = 2^ + ^ +P3) + ^ (fl + <?2 + f 3 ) (1) 

and V, the sum of the three two-particle interaction terms 

V = g (<J(?i - q 2 ) + S(q~2 ~ qs) + S(q 3 - &)) , (2) 

where pj = (h/i)d / dq~j. The potential V represents the idealization of a zero-range 'contact 
interaction'. Because it is comprised of the sum of three one-dimensional delta functions 
and it is self-adjoint, there will be no need for regularization, renormalization or other 
techniques to tame singularities and divergences the full Hamiltonian in the truncated shell 
model (although such techniques can improve the rate of convergence). 

As is usual, we change to scaled, unitless coordinates using the length scale a = y^h/mui: 

q~i = qio and pi = Pih/a. (3) 

Then defining p = (pi,P2,P3) T and q = (q±, q2, <?3) T , the Hamiltonian takes the simpler form 

Ho = ^ (P 2 + q 2 ) (4) 

V = - (%i - q 2 ) + 5{q 2 - q 3 ) + 5{q 3 - g x )) . (5) 

a 

The form Q for Hq suggests the trivial observation that the Hamiltonian of three equal- 
mass harmonic oscillators in one dimension is isomorphic to one isotropic harmonic oscillator 
in three dimensions (forgetting for the moment about any superselection rules). This system 
and its SU(3) symmetry are well-known in the literature. Here we note that the trap 
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Hamiltonian is separable in cartesian, cylindrical and spherical coordinates and therefore we 
can define three classes of energy eigenbasis vectors: 



cartesian: H \rii, n 2 , n 3 ) = fkj(ni + n 2 + n 3 + 3/2)|ni, n 2 , n 3 ) 
cylindrical: H Q \n p , m, n 3 } = hw(2n p + \m\ + n 3 + 3/2)\n p , m, n 3 } 

spherical: H \n r , £, £ 3 ) = Hu(2n r + £ + 3/2)\n r ,£,£ 3 ). (6) 

In these equations, all the n's and £ can be any non-negative integers and m and i 3 can be 
any integer. We define the total excitation number N as 

N = rii + n 2 + n 3 = 2n p + \m\ + n 3 = 2n r + £ (7) 

In the absence of any superselection rules, the eigenspace associated to energy = hu(N + 
3/2) has a dimension of = 1/2(N + 1)(N + 2). The three classes of bases in ^ are just 
different ways of diagonalizing the degeneracies of Hn, the Em subspace of "H. From this 
we can also calculate the dimension of the truncated Hilbert space Y2iLo % w ith maximum 
excitation number N as 1/6(JV + 1)(JV + 2)(N + 3). 

For each of the three types of separable bases in (|6]), the SU(3) symmetry also provides an 
equivalence class of energy eigenbases. Take U(u) as a unitary representation of u e SU(3), 
then any bases of the form 

cartesian: \rii,n 2 ,n 3 ) u = U(u)\rii,n 2 ,n 3 ) 
cylindrical: \n p ,m,n 3 } u = U(u)\n p ,m,n 3 } 

spherical: \v r ,£,£ 3 ) u = U{u) \n r , £, £ 3 ) (8) 

are also complete energy eigenbases for the total Hilbert space %. Further, U(u) only con- 
nects basis vectors with the same total excitation N (note that U(u) may be reducible within 
the degenerate energy eigenspace Hn)- As an example, in the cartesian basis |n l5 n 2 , n 3 ) M , 
the quantum numbers Ui no longer represent particle excitations, but excitations of some 
state corresponding to a particular structure of particle correlations characterized by the 
particular element u. 

Consider the specific element J E SO (3) C SU(3) corresponding to the rotation in config- 
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uration space from particle coordinates q to one set of Jacobi coordinates Jq = r = (x, y, z 



iT. 



J 



/-i- _J, ^ 

1 1 2_ 

l i i 

\ V3 / 



(9) 



Coordinate 2 = l/x/3(gi + ^2 + 93) is the (normalized) center-of-mass of three particles and 
coordinates x = 1/ V2(qi — Q2) and y = 1/ -\/6(gi + 52) — y 2/3g3 are one possible arrangement 
of the standard (normalized) Jacobi coordinates |38j . In these coordinates, the Hamiltonian 
becomes 

#o = y ( k2 + r2 ) ( 10 ) 

V = -§ r (5(x)+5 + f y) + 5 (-ix - gy)) , (11) 

where k = Jp. See Figure[T]for a graphical representation of the Jacobi rotation J. Note that 
since V does not depend on z, the full Hamiltonian still separates in terms of center-of-mass 
and relative coordinates. This means that instead of needing l/6(A^ + l)(A r +2)(A^-|-3) states 
to calculate the ground state energy in a truncated Hilbert space with maximum energy En, 
we only need l/2(N+l)(N+2) because there are only non-zero matrix elements of V between 
vectors with the same quantum number for center-of-mass excitations. 

We will denote the cylindrical H energy eigenbasis in the Jacobi coordinate system as 
|z/, fi, i]} = \n p , m, 71,3} j. In the next section, we will show that these states are a convenient 
basis for irreducible representations of the Cq v symmetry of the full Hamiltonian H, but here 
we note a few additional properties. The center-of-mass motion separates from the relative 
motion so |z/, /1, rj} = \u, fi} ® The wave functions realizing the center-of-mass state \t]} 
in the ^-coordinate are the normal one-dimensional harmonic oscillator states: 



{z\ V } = V>,(*) = vr- 1 / 4 (2V)- 1/2 H v (z) exp (-U^j . 



(12) 



On the plane perpendicular to z that describes the relative motion of the three particles, 
we can define the "hyperspherical coordinates" p and if from x and y using the standard 
conversion to cylindrical coordinates 



P 



\/x 2 + y 2 = Jql + ql + ql- qi q 2 - q 2 q 3 - q 3qi and 



tany? = y/x = Ql J~ 92 ^ . (13) 
V3gi - V3g 2 
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FIG. 1: (Color online) This figure demonstrates that the transformation to Jacobi coordinates can 
be represented as a rotation in configuration space. The % positive particle axes are solid arrows; 
The Jacobi cartesian positive axes are dashed arrows. The z direction points into the middle of 
the solid angle bounded by the all-positive octant and the x and y are orthogonal directions in 
the relative plane, depicted as a gray disk for z = 0, i.e. the center-of-mass at the bottom of the 
potential well. This Jacobi rotation takes place around the yellow (light gray, unlabeled) vector 
and is about 1.22 radians. 

Then the realization of the relative state \v, fi\ as wave functions in the Jacobi cylindrical 
coordinate p and <p are 

Note that only ji = states are not zero at p = and therefore have some probability 
density for configurations in which all three particles are located in the same place. See 
Figure [2] for further explanation of the correspondence between particle configurations and 
Jacobi cylindrical coordinates. 
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V * 



+y / 



1'* 



•PA 



V N 



FIG. 2: (Color online) This figure depicts the relative plane in Jacobi coordinates. The magenta 
(dashed) lines have several interpretations. First, they are the projections of the particle axes on 
the relative plane (see Figure [I]). For example, the line with 99 = ir/6 corresponds to increasing 
separation of particle 1 (depicted as a red circle) from particles 2 (green square) and 3 (blue trian- 
gle), with particle 1 on positive side and the other two co-located on the negative side. Reflections 
across the magenta (dashed) lines also represent two-particle exchanges, e.g. the reflection through 
the line cp = 7r/6 corresponds to the 023 2-cycle. Rotations of ±2-7r/3 are 3-cycles, i.e. three-particle 
exchanges. A rotation of it corresponds to a reversal of parity tt. Reflections across the cyan 
(dotted) lines are a combination of a two-particle exchange and a parity inversion, e.g. 7r<5"23. 

III. SYMMETRIES AND THE SUPERSELECTION RULE 



Expressed in Jacobi cylindrical coordinates, the interaction V takes takes the form [39] 



V 



9 



7T 



(15) 



To understand this functional form, note that the zeros of the delta-functions define three 
planes in full configuration space that intersect the relative xy-plane in three lines (the 
dashed magenta lines in Figure [2]). The locus of points on these planes correspond to 
configurations with two particles in the same position. The rays along (f = tt/6, 5tt/6, and 
37r/2 correspond to the configurations where one particle has a positive g-coordinate and 
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the other two particles have the same negative g-coordinate. The rays with angles <p = tt/2, 
7n/6, and ll7r/6 are configurations with the reverse orientation. 

From this form we can see that the potential V manifests the two-dimensional point 
symmetry that is usually denoted Cq v , i.e. the symmetries of the regular hexagon. The 
finite group C$ v has twelve elements: the identity, two rotations by tt/3, two rotations 
by 27r/3, one rotation by 7r and six reflections of two different types. The group Cq v is 
isomorphic to S3 x Z 2 , where S3 is the permutation group of three objects and would be 
geometrically realized the Jacobi relative plane by C3„. The Z 2 symmetry corresponds to 
the parity invariance of the relative configuration space, which is respected by both H and 
^|40j. We list the group elements and how they act on the relative plane and transform the 
Jacobi coordinates in Table IIHI 

All Cq v group elements leave the hyperradial coordinate p unchanged, but they transform 
if —> ip'. This connection can be inferred from the the geometry of Figure [2j or more 
systematically by using the matrix J to transform the relevant operator. For example, in 
the cartesian particle basis, the three-particle exchange 1 — > 2, 2 — > 3 and 3 — > 1, denoted 
023i is represented in particle coordinate q-space as 



^0 10^ 



<5"231 & 



231 



1 

1 



(16) 



The representation in terms of Jacobi coordinates r-space can be found by transformation 
the matrix with J: 



( _i _vs n ^ 

2 2 w 



J 5*231 J T 



2 2 w 



V 



(17) 



1 J 

which can be recognized as a rotation around the center-of-mass indirection by — 27r/3. 

For convenience, the irreducible representations of Cq v and the characters of the different 
classes of group elements are summarized in Table |TT} These are useful for two reasons. 
First, there will only be non-zero matrix elements of V between vectors in the same type of 
irreducible representation. This can be used to reduce the number of basis vectors that must 
be used to get results accurate up to a certain energy truncation N max . This simplification 



holds for the particular potential (15), but also holds for any two-particle interactions that 
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TABLE I: The first column is the symmetry transformation designated by the corresponding el- 
ement of the point symmetry group of the regular hexagon permutation group Cq v . The second 
column is the same transformation expressed as the corresponding element of S3 x Z 2 . We use the 
notation for S3 permutation group elements such that e is the identity, the 2-cycle a™ switches 
particles i and j, and 3-cycle &ijk switches 1 to i, 2 to j and 3 to The element ir is parity in- 
version. The third column is the rotation or reflection matrix in Jacobi cartesian coordinates, and 
the fourth column is the equivalent transformation of the cylindrical Jacobi coordinate 05 Compare 
with Figure [2] 

depend only on the two-particle separation distances and have parity symmetry, for example 
Gaussian interactions or harmonic interactions. Second, we can use these representation 
to implement the superselection rules. For example, since the cx^'s represent two-particle 
exchanges, the table shows that there are two one-dimensional representations that are 
bosonic under any exchange of particles, A\ and _B 2 , and two one-dimensional representations 
that are fermionic under any exchange of particles, A 2 and B\. The distinction between 
the v4-type and 5-type representations is the sign of the representation of parity inversion 

7T -> C 2 - 



The Jacobi relative cylindrical basis vectors |z/, /z} realized by the functions "i/v^p, (p) (14) 
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-1 


Ei 


2 -2 


-1 
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E 2 


2 2 


-1 


-1 









TABLE II: This the character table for the group Q§ v which can be found in any reference on point 
symmetry groups and we include here for convenience. The entries in the table are the character 
of the representation of the corresponding group element in each irreducible representation, i.e. 
the trace of the matrix that represents the transformation in any basis of the irreducible repre- 
sentation. For example, the trace of the representation of the identity gives the dimension of the 
representation. 

are not actually the basis vectors that transform irreducibly unless \x = (which transforms 
under the trivial identity representation). However, we can define the basis vectors 

1 

1 



|f,A*,+} 



y/2: 



0,/4 + W, -m» 

0,/i} - \u, -fl}) 



(18) 



for any /i > 0. These are realized by the wave functions 

tfUftV) = 1/ ^ „ P^(p 2 )exp(-p 2 /2)cos(/iv?) 



1> v Jj>M 



tt(u + fj,y. 

tt(u + //)! 



p^(p 2 )exp(-p 2 /2)sin(^) 



(19) 



Using ( 19) and the transformations given in Table III, we see the vectors \v, 0}, \v, //, +} and 
\v, fi, — } form useful bases for the irreducible representations of C 6 „. See Table III for the 
exact correspondences. For example, the basis vectors 1^,0} are elements of an A\ bosonic, 
positive parity representation and these vectors are realized by cylindrically symmetric wave 
functions. The \u, 6j, +} vectors (where j is a non-negative integer) also transform under 
A±, so they will mix with the \u, 0} states under the interaction V. Vectors like \v, Qj,+} 
have wave functions with maximal magnitude at three-atom configurations (e.g. = 0, see 
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n ± 
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N.A. 
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BBB, BBX, XYZ 


1 + 


Ei 


Bi 


FFX, XYZ 


1 - 


Ei 


B 2 


BBX, XYZ 


2 + 


E 2 


Ai 


BBX, XYZ 


2 - 


E 2 


A 2 


FFX, XYZ 


3 + 


Bi 


Bi 


FFF, FFX, XYZ 


3 - 


B 2 


B 2 


BBB, BBX, XYZ 


4 + 


E 2 


M 


BBX, XYZ 


4 - 


E 2 


A 2 


FFX, XYZ 


5 + 


Ei 


Bi 


FFX, XYZ 


5 - 


Ei 


B 2 


BBX, XYZ 


6 + 


Ai 


Ai 


BBB, BBX, XYZ 


6 - 


A 2 


A 2 


FFF, FFX, XYZ 



TABLE III: This table identifies the energy eigenbasis vectors \v, 0} for fi = (defined in Section 



2) and the energy eigenbasis vectors for /x > defined in (18) with their corresponding 

symmetry representations and superselection rules. The pattern repeats for vectors with \x > 6. 
BBB (FFF) means three identical bosons (fermions); BBX (FFX) two identical bosons (fermions) 
and one other particle; XYZ three distinguishable particles. 



Figure [2|) and at atom-dimer configurations (e.g. = tt/6). In contrast, any three-fermion 
vectors in A 2 and Bi are realized by wave functions that have nodal lines in the relative 
plane at the atom-dimer configurations (and at the trimer configuration point at the origin). 

The two two-dimensional representations E\ and E 2 do not correspond to three indis- 
tinguishable particles. As with all six representations, they could represent three distin- 
guishable particles, but they also can represent two indistinguishable particles and a third 
distinguishable particle. Assume labels 1 and 2 refer to the two indistinguishable particles 
and must be symmetrized or antisymmetrized. The elements {E, a v , ad, C 2 } form a subgroup 
of Cq v . This point symmetry group C2t, ~ S2 x Z 2 consists of exchanges of only particles 1 
and 2 and parity inversion and therefore can be used to classify superselection sectors for 



BBX and FFX configurations. See Table IV for the character table of this group. The E\ 
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c 2v 
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°"i> 
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1 1 


1 


1 


A 2 


1 1 


-1 


-1 


Bi 


1 -1 


-1 


1 


B 2 


1 -1 


1 


1 



TABLE IV: This the character table for the four, one-dimensional irreducible representations of 
the group C 2v . 



and E 2 irreducible representations of C 6 „ are reducible with respect to the subgroup C 2v , 



leading to the further identifications in Table III 



We should note again that we are assuming that there is no entanglement between the 
spatial wave function and any internal degrees of freedom like spin or hyperfine energy level. 
Additionally, we are assuming that there is not entanglement among the internal degrees of 
freedom themselves. These assumptions imply that these three equal-mass particles all have 
well-defined, separable internal states. When there are three different internal states, we 
call the configuration XYZ. If two internal states are populated, or it is a mixture of bosons 
and fermions, it is either BBX or FFX. Relaxing these assumptions about separability, the 
possible state space for various particle configurations become broader. For example, one 
could place three fermions into superposition of three distinguishable internal states that 
is antisymmetric under pairwise exchange. Then those three fermions could have a spatial 
wave function that was totally symmetric. The possibilities are further extended when 
entanglement between internal state and spatial wave function is considered. 



IV. RESULTS IN TRUNCATED HILBERT SPACE 

In order to find the approximate eigenstates of the full Hamiltonian H = H + V, we 
calculate the matrix elements of the interaction V in the basis states of Hq. The cylindrical 
Jacobi basis states defined in the last section exploit the symmetries of Hq, V and H so 
that, for a given energy level truncation, the fewest number of matrix elements need to 
be calculated. Of course a large improvement occurs simply because the full Hamiltonian 
retains its center-of-mass separability. For a given maximum N, center-of-mass/relative 
separability reduces the order of the model system from 1/6(N + 3)(N + 2)(N + 1) ~ 0(N 3 ) 
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to 1/2(JV + 2)(N + 1) ~ 0(N 2 ) because the effective one-body center-of-mass problem is 
solvable. Additionally, the decomposition of the relative space into irreducible blocks means 
that matrix elements between vectors in different representations are zero, further reducing 
the size of the computation. Although this representation reduction does not change the 
order of the computation from 0(N 2 ), it does reduce the size of computation needed for a 
given truncation size N. For example, for three identical particles with either parity, only 
one out of every twelve states contribute to the calculation. To give a flavor of this, Table 
|V| shows the reduction of the dimension of the truncated Hilbert space for three identical 
bosons and fermions forA^ = 0toA r = 12. 

We calculate the matrix elements in the cylindrical Jacobi basis using their realizations 



in relative coordinate space for p = with the wave function ijj U fi(p } ip) defined in (14) and 



for p > with ip„ }tJ- {p,<p) and (p,tp) defined in (19). The coordinate representations of 
the wave functions are separable in p and ip, and we exploit that to make the following 
definitions: 

{v',0\V\u,0} = -^R(u',0,u,0)A ++ (0,0) 
V 2(7 

{i/,0\V\v,(i,+} = -R(v',Q,v,p)A ++ (0,p) 
a 

{u',p',±\V\u,p,±} = ^-R(u',p',u,p)A ±± (p',p). (20) 

a 



The first two equations in (20) are only relevant for the A\ representation of Cq v associated 



with superselection rules for the cases BBB, BBX or XYZ. The changing factors of y/2 are 



due to the normalization for the p > states defined in (19). 

The functions A ±:iz (p',p) are symmetric in the (p',p) argument and can be calculated 
directly from the integral 

6 r 2n 



, , , v-^ f cos(uo?) cos(uo?) 

A±±(p>, p) = J2 *P ^ - ^) • • (21) 

i=i Jo { sm(// y?) J [ sm(pcp) J 

where each + corresponds to a cosine terms and each — corresponds to a sine term. See 



Table VI for a summary of the results of this integral. Explicit calculation shows that 
74 ±± (/i',/i) will only be non-zero when the particular combinations of {p,±} and {p, ±}' 
come from the same irreducible representations of C6i,. Additionally, we see that states in 
representations A2 and B\ do not feel the two-body interaction potential at all, as expected 
for identical fermions. Wave functions in these sectors all have nodal lines corresponding to 
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TABLE V: The total excitation number N identifies the energy eigenspace with En = huj(N +3/2). 

The dimension of the truncated Hilbert space up to and including that energy level is Dn = 

1/6(N + 1)(N + 2)(N + 3). We define as the dimension of the truncated Hilbert space with 

maximum excitation N but constrained such that center-of-mass motion is in the ground state. 

In other words, for the Jacobi cylindrical basis, D N is the dimension of all relative states of 

distinguishable particles such that r/ = and 2v + < N. The other dimensions are: Df+ for 

three identical bosons with positive parity, for three identical bosons with negative parity, 

and similar definition with fermions for D 3 - F+ and D 3 ~ F_ . 

N N 

configurations where two particles overlap, and are therefore antisymmetric under reflection 
across this line [41J. Therefore, the three indistinguishable fermionic relative states of the 
form \u, 6j + 3, +} and \u, 6j, — } are eigenstates of both H and H. 

To analyze states in the other four irreducible representations, we must calculate the 
radial integral: 

W^M) = ,. V Z', = / rfpp^'+WKV) e- p2 - (22) 
iiJiu + u)\{i>' + u'y. Jo 
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FFX+ 


A--(6j + 3±l,6A; + 3±l) 


= 3(- 





TABLE VI: This table summarized the results of the angular integral (21) for the different ir- 
reducible representations of C% v and C 2v . The letters of the third column describe the 'most 
indistinguishable' three particle configuration with that representation and the corresponding par- 
ity. For example, "BBB+" in the first line can also correspond to BBX and XYZ states with 
positive parity. 



Using the substitution p 2 = t, this can be put into the standard form 



which can be evaluated as 



27iV(i/ + + //')! Jo 



dtt^^L* (t)L» f (t) e 



(u' + //)! T(u + i+i^)r(i+^) 



27rV z/!z/ / !(z/ + /i)! 
x 3 F 2 (-1/', i±^/,t 



r(i±^; 



^;-u+ i=^,/i + 1;1) 



(23) 



(24) 



The function 3 F 2 ( 



is a generalized hypergeometric function. The lack of symmetry 



between primed and unprimed coordinates in (25) is deceiving. One can either use the 



definition (22) or the properties of ^F 2 { 



to show that R(v', //, v, \i) = i?(i/, p, z/, fjf) 



Combining the results of Table VI and equation (|25j) to calculate the matrix elements 
(20), we can proceed with some numerical results. Figures [3] and [1] depict the relative energy 
levels calculated up to N max = 2v + p = 30 (i.e. E™ l ax = 31hu) for the two indistinguishable 
boson sectors. This corresponds to a 51-state basis for the positive-parity bosonic sector Ai 
and 40 state basis for the bosonic B 2 sector. Note that in these energy level diagrams, all 
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FIG. 3: (Color online) Solid black lines are the 40 energy levels with symmetry corresponding 
to the totally symmetric, positive parity A\ representation of C§ v for N max = 30. Energies are 
measured in units of Hlj. The dashed red lines are the exact result for universal states in the 
limit g/a — > oo and correspond to the energy levels and degeneracies of the non- interacting FFF+ 
A2 representation. This representation contains the ground state of the BBB, BBX and XY Z 
configurations that is a trimer in the case of attractive interactions. See text for more details. 

apparent crossings are actually avoided crossings. Incorporating center-of-mass motion only 
adds a shift of fko{r] + 1/2) to each of these levels. 

Looking at Figure [3j three different kinds of states are identifiable. First, the states whose 
energies are are approximately horizontal for g/a ~ E^ l ax are three-atom states. This can 
be most easily seen for the low energy states, where as \g/cr\ becomes much larger than 
hw the bosonic states become effectively fermionic, i.e. they approach the exact energy and 
degeneracy expected for the non-interacting, postive-parity three-fermion Bi sector. The 
wave functions for these states all have maximum amplitudes in the sectors of relative space 
corresponding to three disparate particles, specifically ip = 0, ±vr/3, and ±27r/3. In other 
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FIG. 4: (Color online) Solid black lines are the 40 energy levels with symmetry corresponding 
to the totally symmetric, negative parity B2 representation of C§ v for N max = 30. Energies are 
measured in units of hu. The dashed red lines are the exact result for universal states in the 
limit g/a — > and correspond to the energy levels and degeneracies of the non-interacting FFF- B\ 
representation. 

words, in these spectra we can observe the Tonks-Girardeau phenomenon of fermionization 
of bosons with strong interactions . These states can be considered universal in the 

unitary limit of \g/cr\ ^> hu because the energy levels only depend on the trap parameters 
and not the interaction [HI [26] . For higher energy states near the cutoff, the numerical 
agreement between the bosonic state and the corresponding fermionic limit becomes less 
precise, although the bands still have the correct degeneracy. This systematic effect is due 
to the truncation of the Hilbert space, and occurs even though the interaction potential is 
renormalizable in one-dimension. Methods from nuclear physics exist for transforming the 
Hamiltonian into an effective Hamiltonian to account for this effect, e.g. the Block-Horowitz 
approach used in [TO] or the Lee-Suzuki method from [29], and these will be explored in 
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subsequent work. 

For g < 0, the three-atom 'gaseous' states would not be stable in Figure |3j but could 
decay into 'liquidlike/solidlike' states [2j consisting of the two other types of states: trimer 
states and atom-dimer states. The trimer state only occurs in the A\ sector because only 
/i = states can have non-zero probabilities amplitudes for all three particles to be in the 
same location. Above this state are the atom-dimer states, which diverge as g decreases from 
zero and the dimer becomes more tightly bound. For \g/cr\ < E^ l ax , the energy eigenvalue 
diverges roughly quadratically, and the difference between the different atom-dimer states 
is just the energy cost of exciting the atom and/or dimer to a greater value of v. For 
\gjo\ ^ E^ l ax , the atom-dimer and trimer states appear to diverge linearly, but this is again 
a limitation of the truncated Hilbert space. As one increases the number of basis vectors in 
this section, the expected quadratic divergence holds for larger values of of the interaction 
parameter g. Additionally, more atom-dimer states appear and the spacing between those 
levels for a given g < becomes more exactly 2hu. 

A couple other interesting features are the stepped cascades of gaseous states for increas- 
ing negative values of g and the divergent states with g > 0. The divergent states for g > 
are again an artifact of the truncation. As the cutoff is raised and additional energy lev- 
els are added, the lowest divergent states bend down to the fermionic limit. On the other 
hand, the cascade structure is not an artifact of the truncation. As higher energy levels 
are added, the steps of the cascading energy level become more regular for low negative g. 
These structures depends on coupling between /i states of different relative (configuration 
space) angular momenta and are absent when that coupling is removed. 

As a final comparison, we also include Figure [5] which depicts the spectrum for the E\jB\ 
representation, corresponding to an FFX configuration, such as two spin-up fermions and 
one spin-down fermion. 

V. CONCLUDING REMARKS 

Some extensions to more particles of the method described here are straightforward. 
For example, the representation theory of permutation groups in three dimensions is well- 
known, and an extension to four particles in one-dimension can be accomplished using the 
octahedral symmetry of the relative configuration space. Further extensions to more than 
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FIG. 5: (Color online) Solid black lines are the 80 energy levels with symmetry corresponding 
to the two-component fermion with negative parity in the E1/B2 representation of Cq v /C2v for 
N m ax = 30. Energies are measured in units of Hlj. This is the ground state sector for the two- 
component fermion. The dashed red lines are the exact result for universal states for both the 
FFF+ A2 and the FFF- B\ representations, showing the emergence of energy levels corresponding 
to 'anomalous' parity states in the unitary limit. 

five particles in one dimension require using the less-familiar (at least from molecular and 
solid state physics) point groups in higher dimensions, which have been classified by Coxeter 
and others. Since the three-body problem in one-dimension carries so many properties of 
the three-dimensional system, such as universality and bound state structure, this line of 
investigation seems promising. 

When extending this method to higher particle numbers and higher dimensions, for 
example N particles in d dimensions, one difficulty is that the realization of the permutation 
group Sat in 0((N — l)d) no longer commutes with the 0((N — 1) ^-representation of the 
0(d) symmetry of the two-particle interaction, a symmetry that is also no longer a discrete 
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symmetry. Contrast this with the present case: the 0(2) realization of S3 in the Jacobi 
plane commutes with the 0(2) parity representation of the discrete Z 2 = 0(1) symmetry 
of the interaction. However, we do not find this task hopeless, and note that even in d- 
dimensions, there is always 0((N — l)d) representations of particle permutations. Such 
geometric realizations of the symmetries of the relative configurations of three particles 
appear in the construction of Talmi-Moshinksy brackets used in [151 HH] and the perturbative 
expansion solution of Faddeev equations in [6] . 

Additionally, in higher dimensions there are the difficulties of renormalization of the 
zero-range interaction, but methods from effective field theory and/or regularization have 
previously been used to study these interactions in truncated Hilbert spaces and their re- 
sults should be extensible. Similar approaches will also be necessary if genuine multi-body 
interactions are included; for example, a term proportional to 5(p) in the three-body, one- 
dimensional case. 

For anisotropic traps, the rotational symmetry is gone, but there can still be reflection 
symmetries whose representations can be exploited for constructing symmetry adapted bases 
using these methods, as shown by [2T]. Similarly, even with particle of uneven masses, 
discrete symmetries can help reduce the problem in a numerically more efficient way. 

As a final comment, the same symmetry techniques that provide efficiency through sepa- 
rability and reducibility in the present work should also be useful for calculating the entangle- 
ment among atoms using truncated Hilbert spaces. The coefficients that connect the Jacobi 
cylindrical basis can be calculated using a variety of symmetry-based methods. In the case 
of indistinguishable bosons, these coefficients can be used to find the distance from the en- 
ergy eigenstate to a separable state as quantified by any atom-dimer bipartition. In the case 
of indistinguishable fermions, the entanglement can be compared to the least entanglement 
one can expect in an antisymmetrized state. In either case, we believe that entanglement 
spectroscopy of trapped few-body systems could be an interesting testing ground for the 
interplay of particle superselection roles and the emergence of composite systems. 
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